RADIO SCIENCE lournal of Research NBS/USNC-URSI 
Vol. 69D, No. 7, luly 1965 

Calculated Curves for Groundwave Propagation 

Over Inhomogeneous Earth With Pronounced 

Topographical Features 

K. Furutsu 
National Bureau of Standards, Boulder, Colo. 

(Received December 21, 1964; revised January 28, 1965) 

Several numerical examples are presented for radio propagation over a three-section path. These 
are based on formulas which were derived earlier [Furutsu, 1957a, 1957b, 1959, 1963; Furutsu, Wilker- 
son, and Hartmann, 1964] in which the height and the electrical properties of each section may be 
different. First the path is assumed to be a land-sea-land path and a sea-land-sea path with no height 
difference; the attenuation using this model is obtained both for the flat earth and the spherical earth 
and the values are compared to show the asymptotic approach of these two values at short distances. 
Second, the effect of the height and width of a mesa is illustrated and the attenuation values are com- 
pared to the values obtained when the obstacle is a knife-edge of the same height. Next the effect 
of varying the receiver height is presented when the path is either a sea-land-sea path or a land-sca-Iand 
path. Finally, examples of the attenuation caused by both a channel and an island having bluffs are 
given which illustrate the recovery effect as well as the diffraction loss in the vicinity of the coastlines. 



1. Introduction 

This paper presents some numerical results obtained assuming a three-section inhomogeneous 
earth which may be either spherical or flat. It is a continuation of a paper which assumed a two- 
section earth [Furutsu, Wilkerson, and Hartmann, 1964] . These two papers are based on a theory 
of radio wave propagation which was developed previously [Furutsu, 1955, 1957a, 1957b, 1959, 
1963]. The basic terrain consists of three sections of different heights and electrical properties 
with a ridge at each section boundary. This terrain is illustrated in figure 1. Using this model 
the field may be found over an inhomogeneous earth path consisting of three sections, such as a 
land-sea-land path having a bluff at each coastline, a two ridge path, or a path including a plateau, 
etc. 

The basic formula of field strength can be given in a form of mukiple residue series, which 
reduces to the ordinary Van der Pol and Bremmer series in the special case of a homogeneous 
earth. Also the corresponding flat earth formula has been obtained in an analytical form for the 
limited case where the earth's surface is smooth and both the transmitter and receiver are on the 
ground [Furutsu, 1955]. This flat earth formula may be used for short distance paths when the 
spherical earth formula is poorly convergent. 

Figure 1 shows the basic terrain used in this paper, which consists of three sections of different 
earth radii, a2, a4, and ae, and different propagation constants, A:2, A:4, and /ce, respectively, with 
ridges of radial heights, as and a^ (measured from the earth's center to the tops of the ridges), at 
two boundaries of the three sections. The transmitter and the receiver are located at the points 




Figure 1. The geometry of the basic terrain for equation 
(L2). 
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Xi and X7 whose radial heights are Zi and Z7, respectively. Also, the path lengths along these three 
sections are respectively r2, r4, and re, which are assumed to be measured along a mean earth 
surface of radius a. 

The attenuation coefficient A is defined in such a way that, if E is the field strength to be 
obtained, 

E = 2AEo, ^, = __— i— _-e-Mr.+-4+r6) k, = 2TTlK (1.1) 

47r(r2 + r4 + re) ^ ^ 

where X is the wave length in free space. Thus, j^o may be regarded as the field strength in free 
space excited by a dipole of the moment (47r)~^ (in Gaussian units). 

Then, according to formula (38) in [Furutsu, Wilkerson, and Hartmann, 1964], the attenua- 
tion A in this particular case becomes as 



A = S,„ u, t, V77(C2 + C4 + C6) e--/y,,(ri2) (^2 " gD'^e"^"^- {q^;Sys2)ftSys4) - q^i{TM)fdy22)} 

X (72 - 74 + ^2 - ur\t^ - d)-'e-^^^^^ {94//,(y54)A(r56) - q6f;,{y5e)My^d} 

X (74 - ye + ^4 - ^6)-H^6 - ql)-'e-i^''%{y76), (1.2) 
where C2, C4, and Ce are the numerical distances corresponding to the distances r2, r4, and re, i.e., 

C2 = (r2/a) ik^a/2yl\ c, = (rja) {k,al2yl\ Ce - (rela) {hal2yl\ (1.3) 

and 

yij=yi-yh (1.4) 

where ji (i= 1, 2, . . .,7) are the numerical heights defined by 

n = Hz, - a) {2lk,ay'\ 77 = h{z, - a) {2lk,ayl\ 

yi = h{ai-a){2lk,ayl\ i = 2, 3, . . ., 6. (1.5) 

Here the set of ^m im = 2, 4, 6) stands for the roots of the equation 

W'{tm)-qmmtm) = 0, (1.6) 

where W'it) is the first derivative of the function W{t) defined by 

r(- t) = (77t/3)l/2e-^-2-/3//(2) g ^3/2^ 

and 

qm=^-i(k^al2yl\hlk'J, (1.8a) 

with 



(1.7) 



ikml\^k^ — kl, Vertical Polarization 
(1.8b) 
A;2/\/A;^ — A:J, Horizontal Polarization. 
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Also 

ftjy) = mtm-y)mt,n), 

flm 'y) = ^'(tm - y)/W '«„,) = g- ' W'(t,n - yWiU), (1.9) 

where ftjiy) is the ordinary height gain function. From these definitions it follows that 

/'n.(0)=/;jO)=l- (1.10) 

The corresponding attenuation coefficient in the flat earth approximation can be given in terms 
of the following functions : 

g{z, )8, a) = f^)'e('+^^('+«^)l-^f"d<,e-'? \" dhe'^i t^ e-'ldh, 

Xiz, /3, a) = (-^Y e('+«^+''^'-^f " dt,e-'^ f " e-'ldt-, f " e->Mh, (1.11a) 
Wtt/ J, Jut, Joct, 

g^^\z) = [dldzYS'(z), 

<?'»)(z, a) = idldzYS{z, a) = 2(l + a2){2^<"-"(z, a) + (n- l)#<"-2»(z, a)} -(2/\/^)a»-'r$'(''-i)(a2), 

X<'>(z, )8, a) = 0/a2)x(2, )8, a) = 2(l + a2 + y82)zx(2, i8, a)~&\fii-)8{ca)8(fiz). (1.11b) 

with the relations 

g{z,a) = S' iz) S'iaz) - ,^ (oz, «-•), a > 

(^(z)^'(/3z, a)- S'iz, ^, a)=i(aPz)S'{Pz, /3->)-^ (a/3z, a"', /3-') 

= X()8z,/3-', a), a,/3>0. (1.11c) 

The attenuation coefficient, A, is given in appendix A for the general case. In the special case of 
K ~ Ky which may be practically most important (land-sea-land and sea-land-sea paths), the attenua- 
tion is given by 

— S {iVm2p4, V(m4 + me)lm2) — ^{iVmep^, V(m2 + m4)lm^ 

-f S {iy/m2p4y V/^6/^2, ^/mjme) + S^{ivrneP4, Vm2/m6, \/m^lm<^ 



- {k'Jk'2) {{g {iVp4{m4 + /Me)) — ^' (iVm6P4, VW^)) S{i^rn2P2 



-\-{8{i vp4ijn4 -f 7712)) — ^^ (i V/n2P4, vnHlm^)) ^ {i Vm6P2)} 
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-\-{S{i vp2{m4 + me)) — S{i \rrup2'> Vmelrru)) 6^^\i ym^pi) 



-^{8ii Vp2(m4 + 7M2)) — ^ (i Vm4p2, \m2lmA)) 8 ^^\i ymep2) 

+ 2(x(i Vm4p2, vmolmA, \/m2lm^ + i V/714P2 X^^H^ ^rn^p^, wmolrrU', \Tn2lm4)) 

— 2i(V(/7i4 + m6)p2 ^^^\i yp2{m4 + me), wm^lim^ + me)) 

4- V(m4 + m2)p2 (^^^^\i Vp2(m4 + m2), Vme/(m4-l-m2)))j 

r / A:' \^ A:' — A:' If 

+ |2(^r^) + ^^+1^ ~ "^ J I ^'^^ ^'^2(^4 + me) , Vm2/(m4 + me)) 

+ (^ (i Vp2(m4 + m2), Vme/(m4 + m2)) J 

7'_L7M ) — tan-^A/ — ^(i Vp^)H-x(^ Vm4p2, Vme/m4, Vm2/m4)[ . 

A:4 H- /Cg/ LTT \ m2m6 J JJJ 

(1.12) 



-2'—^ 



Here 



ni2 = r2lr, m4 = nlr, me = relr, r=r2 + r4 + re. 



(1.13) 



and p2 and p4 are the total Sommerfeld's numerical distances measured by the ground constants 
fe and A:4, respectively. 
Using the identities 



(iVm2p4, vme/m2, V m4/me) — <^ (i V m2p4) ^(iVmep4, vmjme) 



-h S' (ivmep^y Vm2/me, Vm4/m2)— S {iwmQp^} & (iN mip^-, vm^lm^) 



= — X(iVmep4, Vmjme, Vmjme) — x(^ Vm2P4, Vm6/m2, VWm^) 



~X(^^m4P4, Vm6/m4, V m2/m4)— ^(i V m2p4)^(i V m4P4) ^(i V mep4), 



(1.14) 



the formula (1.12) in the special case of k2i=kQ) = k'^ can be shown to reduce to that for a homoge- 
neous earth of the ground constant A:2. The analytic expansions of the functions used are treated 
in appendix B. 

The Sommerfeld numerical distances p2 and p4 introduced in (1.13) and the numerical dis- 
tances C2, C4, and ce defined by (1.3) are related by 



P2 = iqlic2 4- C4 + Ce), p4 = iqlic2 + C4 + Ce). 



(1.15) 
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When the frequency is low enough so that €m is much less than 18 <Tm//, the Sommerfeld nu- 
merical distance is given by 



p„ = (577/27) (fla,„) (r., + r, + r«). 



(1.16) 



In the above, / is the frequency in Mc/s, (Jm is the conductivity in millimho/m and r? is the distance 
in km. 

In concluding this section, it may be noticed that the theory is vahd independent of the incident 
angles of wave on the boundary lines of sections, except in the immediate vicinities of the bound- 
aries. On the other hand, it does not take into account the possible wave components which are 
reflected several times between the two boundaries of sections (fig. I) and are finally propagated 
to the receiver. In the case of mixed paths (smooth earth), the attenuation also has been obtained 
in terms of convolution integrals [Feinberg, 1944, 1946, 1959; Bremmer, 1954; Wait and House- 
holder, 1956, 1957; Wait, 1961], and is usually computed by a numerical method. 



2. Three-Section Mixed Paths 

In figure 2a, the flat earth value (broken line) of the attenuation is compared to the spherical 
earth value (solid line) for a sea-land-sea path (curve 1) and also a land-sea-land path (curve 4). 



<- 10-2 




1000 10 20 50 100 200 500 1000 

P P 

Figures 2a and 2b. Comparison of the flat earth values of the attenuation to the corresponding spheri- 
cal earth values for various combinations of the ground constants, when r-z = r^ == ru and the sea has 
the constant (a) q = 0, (b) q= LSe"'''/^. 
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The graph is plotted against the total Sommerfeld numerical distance, p, measured using the con- 
ductivity of the land, ctl, (€l ^ o-Jf) when ri^r^ — fq. The sea is assumed to be perfectly conduct- 
ing and the ground constant of the land is assumed to be q = 30e~^'^/'^. In this figure, the curves 
2 and 3 are also shown when the paths are reduced to the mixed paths of two sections caused by 
putting qz = Qa. 

In figure 2b, the attenuations are found for the same paths on the same conditions as in figure 
2a except that instead of a perfectly conducting sea the ratio of the sea conductivity, C7s, to that of 
the land, ct/,, is assumed to be 400 (this gives the ground constant ^ = 1.5e~^^/^ for the sea parts). 
As in figure 2a, the curves are also shown for the associated two-section paths. 

As expected the spherical earth value approaches the flat earth value asymptotically at short 
distances in both figures. 

Figures 3a and 3b illustrate a set of curves of the attenuation and phase for land-sea-land paths 
when the earth is flat and cts/cta^ 400. They are displayed versus the total Sommerfeld numerical 
distance measured by using the land conductivity cr/,, for various values of the parameter A^= (r2 + re) 
/r4 in the particular case of rz — r^. 

It may be noticed from figure 3b that, in the case of the flat earth approximation, the phase 
delay of the attenuation coefficient for a land-sea-land path can exceed tt for p > I and A^ <^ 1, when 
the finite conductivity of the sea is taken into account and the corresponding numerical distance of 
the sea part becomes of the order of magnitude of I or greater. A similar situation can also be 
seen for a land-sea path, a mixed path of two sections, on the same conditions. For comparison, 
the attenuation and the phase in the latter case of land-sea paths are also displayed in figures 4a 
and 4b on the same condition, o-s/ct/^ = 400, as in figures 3a and 3b, where N = r2/r4, the ratio of the 
land part to the sea part, and p is the total Sommerfeld numerical distance measured by using the 
land conductivity (refer to the equations with rQ = in fig. 3a). 
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Figures 3a and 3b. Flat earth values of (a) the attenuation and (b) the phase delay for land-sea- land paths 
versus the total Sommerf eld's numerical distance p (measured by the land conductivity ctl) for various values 
ofN = (r2 + r6)/r4, when cr J <Ji^ — 400 and i^/re = 1 . 
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Figures 4a and 4b. Flat earth values of (a) the attenuation and (b) the phase delay for land-sea paths, versus the total 
Sommerf eld's numerical distance p (measured by the land conductivity cti) for various values ofN, the ratio of the land 
part to the sea part, when aJai=4-00. 




Figure 5. The geometry of the mesa for figures 6a and 6b. 



■ RIDGE VALUES 

■ MESA VALUES 



Y34=Y54 = Y76 = 




0.05 0.07 0.1 0.2 0.5 

NUMERICAL DISTANCE OF MESA (C4) 




12 3 4 

NUMERICAL HEIGHT OF MESA 



Fl(;uRE 6a. The attenuation caused by a mesa versus 
the mesa iength compared to the attenuation caused 
by a ridge of the same height. 



Fl(;URE 6b. The attenuation caused by a mesa relative 
to the corresponding ridge value versus its height for 
various values of the mesa length (the absolute values 
are displayed). 
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3. Effect of a Mesa 

In order to illustrate the effect of a mesa the terrain is considered to be as shown in figure 5. 
The distances on each side of the mesa, C2 and Ce, are both equal to 1 and the transmitter and the 
receiver are both on the ground. In figure 6a the attenuation is illustrated versus the numerical 
width of the mesa, C4, for its numerical height, 732 = y56— 1 and 5. 

Also shown in this figure is the value of the attenuation which would be obtained if a ridge 
(broken line) of the same height as the mesa (sohd line) were placed at the midpoint of the path 
with the same total length. It definitely shows the obstacle gain (the higher the obstacle the smaller 
the attenuation), and also that the effect of the mesa width becomes greater with increasing height. 

The latter situation is also illustrated in figure 6b, where the attenuation relative to the cor- 
responding ridge value is shown versus the numerical height of the meas for various values of its 
width with the same values of other parameters as in figure 6a. ^ 



4. Effect of Varying Receiver Height 

In the previous paper [Furutsu, Wilkerson, and Hartmann, 1964] the variation of the field 
strength was illustrated when the receiver was raised and the path was a mixed path of two sections 
(sea-land). A similar situation is illustrated in figure 7 for a mixed path consisting of three sections. 

In figure 8a, a path is selected corresponding to a land-sea-land path having the numerical 
distances C2 = 10, C4=1.0, and C6 = 0.7 with the ground constants q2 = q6 = — iS0 and q4=—i 10"^ 
and, keeping the transmitter (xi) on the ground, the numerical height yye of the receiver (xj) is 
changed; the attenuation is shown by a solid line, while the broken and the chain lines are the 
associated mixed path values of two sections and the homogeneous earth values for ^2=~i50, 
q4 = q6 = ~i 10~^, and q2 = q4=^q6 = ~i 50, respectively. The two optical boundary points of this 
path are found to be at y76 = {0.1)^ and (1.7)^. However, as is seen, the three section values ap- 
proach the two section values and then the homogeneous earth values at much higher points than 
these optical boundary points. It also shows some oscillation, which may be interpreted as an 
interference of the principal wave with the waves induced at the section boundaries. 

Figure 8b shows a similar illustration for a sea-land-sea path with the same values of the 
parameters involved except that q2 = qQ = ~i 10"^ and ^4 — ~ i 50. 



5. Examples of Field Strength When the Obstacle Is an Island or a Channel 

In figure 9, a sea-land-sea path (island) is considered where the height of the island is 500 m. 
The transmitter and the receiver are both assumed to be on the ground, the wave frequency is 1 
Mc/s, and vertical polarization is assumed. The island is 200 km from the transmitter, and 
the electrical constants of the land are eL=^ 4, (Tl= 10 millimho/m and those of the sea are €§ = 80, 
0-5 = 4000 milhmho/m. 



' When the frequency is 300 Mc/s, the value C4 = 0.5 corresponds to the distance 14.2 km and the value Y^z — b to the height 238 m. 




Figure 7. The geometry and the optical boundary points 
A \ for figures 8a and 8b. 
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Figures 8a and 8b. The variation of attenuation versus height over mixed paths of 3 sections along with the connected mixed 

path values of 2 sections and the homogeneous earth values. 

(a) a land-sea-land path, (b) a sea-land-sea path. 




Figure 9. The geometry of an island having bluffs for 
figures 10a and I Ob. 



Figure 10a shows the field strength in dB (0 dB = 1 mV/rn) versus distance (km) for several 
lengths of the island, when the 2|£'o| [defined by (1.1)] is assumed to be given by 300/(r2 + r4 -h re) 
mV/m.^ Figure 10b shows the corresponding phase delay in cycles. Here it may be noticed 
that the phase of £"0 defined by (1.1) is the phase delay corresponding to the distance measured 
along the surface of a mean earth radius a (fig. 1). In figure 10b and also in the following figures 
lib and 13b, the relative phase is displayed when this mean surface is taken at the sea level. 

Figures 10a and 10b show the recovery effect as weU as the diffraction loss in the vicinity of 
the bluff at the second coastline. The effect of the bluffs may also be seen by comparing to figures 
11a and lib, in which no bluff is assumed with the same values for other parameters as in figures 
10a and 10b. 



^ This value is the field strength when the transmitter of a vertical electric dipole and the receiver are both placed on a perfectly conducting plane earth and the 
radiation power is IkW. 
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Figures 10a and 10b. An illustration of (a) the field strength and (b) the phase delay versus distance for island paths. 

(The broken lines are the expected continuations of the solid lines for the 3d section.) 
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Figures 11a and lib. (a) The field strength and (b) the phase delay versus distance when the island in figure 10 has no 

bluffs. 
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Figure 12. 



77?^' geometry of a channel for figures 13a 
and 13b. 
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Figures 13a and 13b. An illustration of (a) the field strength and (b) the phase delay versus distance for channel 

paths. 



In figure 12, a land-sea-land path is illustrated, where the height of the land sides are again 
500 m and the other conditions are assumed to be the same as in the preceding case of an island 
path. The solid lines in figure 13a show the field strength along the ground surface, while the 
broken hne shows the field strength 500 m above the sea section, this being the same height as the 
land. Generally, the latter values would be lower than the former (up to about 6 dB) except for the 
vicinities of the bluff where the diffraction loss becomes predominant, since a smaller contribution 
is expected from the reflected wave from the sea surface. Figure 13b shows the corresponding 
phase delay. 

Finally, it may be noticed that, in figures 10, 11, and 13, the possible contribution of the sky- 
wave by the ionosphere is not taken into account, which becomes predominant at large distances 
[Wait and Walters, 1963]. 

6. Conclusion 

Numerical results based on previous papers [Furutsu, 1957a, 1957b, 1963] are presented for 
propagation over an inhomogeneous earth consisting of three sections with different heights and 
electrical properties; a ridge may also be located at each section boundary. This is illustrated 
in figure 1. 

The flat earth formula [Furutsu, 1955] is used to show the asymptotic approach of the spherical 
earth values and the corresponding flat earth values at short distances, although the flat earth for- 
mula is available only when the earth's surface is smooth and both the transmitter and the receiver 
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are on the ground. A set of curves is illustrated for the flat earth values of the attenuation and the 
phase delay of land-sea-land paths. It shows the fact that, when the finite conductivity of the sea 
is taken into account and the Sommerfeld numerical distance of the sea part becomes of the order 
of magnitude of 1 or greater, the phase delay of the attenuation coefficient can exceed n although 
it can not in the case of a homogeneous flat earth. 

The effect of a mesa, as illustrated in figure 5, is considered and the values of the attenuation 
are compared with the corresponding values obtained when the obstacle is a ridge of the same 
height. The mesa and the knife edge both exhibit the eff'ect of the obstacle gain, and the effect 
of the mesa length is found to become more predominant with increasing height. 

The variation of attenuation is also illustrated as the receiver is raised when the path is a 
land-sea-land path or a sea-land-sea path; this situation is illustrated in figure 7. As may be seen 
from this figure, the field strength computed over a three-section mixed path is expected to depend 
mostly on the sections closer to the transmitter as the receiver exceeds the heights of the optical 
boundary points. Thus the attenuation is compared with the mixed path values of two sections 
and further with the homogeneous earth values. However, the mutual approach of these values 
is found at much higher points than the corresponding optical boundary points. 

Finally, in figures 10a, 10b, 13a, and 13b, the field strength and the phase delay are illustrated 
for a sea-land-sea path and a land-sea-land path, respectively, when the land is 500 m above the 
sea; both the transmitter and the receiver are on the ground and the frequency is 1 Mc/s. The 
recovery effect as well as the diffraction loss in the vicinities of the bluffs at the coastlines are 
observed in these figures. 

Although the present theory is valid only for terrains with sharp vertical faces, as illustrated 
in figure 1, it may be available for other cases, in the same way as the knife edge model has been 
found to be useful for actual ridge diffraction of VHF waves. 



The author thanks R. E. Wilkerson, Perry H. Elder, and Raymond F. Hartmann for their help 
in obtaining the numerical results for this paper. 



7. Appendix A 

Putting 



A = l- V7Tki(r2 + r4 -f re)l2 e'^^'^B, 
Pz =- UhlKfhr^n, p4 = - i{hlk',fkynl2, pe=- i{krlk'^fk,r,l2 

[these notations p^'s are to be distinguished from the same notations used in (1.12)], then 



B = 



/t 1 Kn'^fi 



k'^{k[ -h kl2) {k'^ + A:^) W V m2m% 



tan 1 A / — g {i Vp4lm4) 



— 6 {i v p^m^lrrtA^ V(m4 + m^\mi) — fj {i vpimjTru, V(m4H- m-iMm^) 



+ lo {i Vp4m2/m4, \frn^lm2, vmjme) -\- S H Vp4rnjm4, Vrn^frne, Vm^m^) 



~^ ik' -\- k')lk' -\- k') ^^^^ ^ ^^ Vp4(77i4 + mQ)lm4) — 6 (i VpJnJm, Vmjme)) g {i Vp^) 



+ k2{ & {i \/p4{m4-^m2)lm4) — 8 {i Vpm^Jrru, Vm4lm2)) S'(i Vp^) 
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- K f^a V^) /; a V^) //(i V^)} + (^7^) [§ /; a V^J^) 



K2 



Vp2/m2) I 



+ 






(Al) 



k, 



^f— ^) {A:^ /Ai Vp2(m4 + //i2)/m2) ^''ii Vp^) 



-k'^fiiiVpzm^lmz, V{m4 + m2)/m6)} l^/LV l^I //(i Vp6(/7i4H- /n2)/^6, Vm^limA-^m^)) 



% ^2 



A^4 T A)(j 



M^6^4 "^2 ) /;/• - 



— k^SiiVpemzlnifi, VimT-K^e)/^)} ^A^^;;;^) //(i Vp2{m4 + me)lm2, VW(^i7-K^) 



/:, 



fin IVry 



J , J j {A:^ ^'(i Vp^) /? (/ \^p^mjm, VmJniA) 



7—7) {A:2 '^'(^* Vp^) //(i Vpzmjmz, \Tri2jm4) — k^ ^{i ^ pim^lm^, Vmjme, Vmjm)} • 



8. Appendix B 



The function (ff (z, /3, a) can be given in terms of the x-function using the relation (1.11c). 
Therefore, in this appendix, the analytical expansions will be given only for the functions // (z), 
io(z, a) and xU, «, Z^)- 

When a and fS are real positive, the definitions of these functions in (1.11a) give the following 
expressions: 

77 

X(z, 13, a)==xiz, a, ^) = e^'''' P>{z, a) + e°^^^'^(2, |8) 

+ X(0, a, I3)€^'^'"*i^'^''-e'"'^l''^''i^{z)-x{z, «, /3). (Bl) 

Here 



X(0,a, /3) = -tan-> 



(l-a/3)Vl + a^ + /3^ + a^(Q; + ffl 
(a + /3) VTT^^^:^^ + a)8(a/3 - 1) J ' 






X(z, a, i8) = f-^j e('+«'+^^)^M c^^ie'l e-'y^2 e-^3rf^3. 



(B2) 



767-936 0-65— 9 
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The power series expansions of f) (z), ^ (z, a), and x(2, «, /3) with respect to z are given by 

n=l 

c„ = (re-V2)-'2='c„-i, ci = {2lV^)z, (B3) 

#(z, a)= J C^n, 

dn = n-'{{l-\-a'')z'dn-i-^fn}, (B4) 

00 

X(z, a, j8)=21 •^^' 

n=l 
Xn - (/I + V2)-H(l + ^2 + ^ V^^-1 + yn} , 

(B5) 

an = {n — V2)-iQ:Van-i, bn = {n— V2)~^j8 V6n-1, 
ai = 6i = (2/V^>)82^ 

While (Bl) gives the power series expansions with respect to z, the asymptotic expansions 
for large values of z are found to be given by the following equations: 

n— 1/2 1 

g{z)-^gn. \z\»\, |arg(z)|<37r/4, ^„- ^ gn-u So = -j^^^ (B6) 

1 » 

fAz, a) ~ y yn(a/ V 1 + a\ R), 7? ^ vTT^ z, |/?| » 1, | arg («) | < 377/4^ 

V7!-(l + a2)z^^^ 

X{z, a, 13) ~ -^ ^- f h„, 



R = Vl + o^T^ z, \R\»l, |arg («)| <37r/4, (B7) 
with 

ho=i{az)^'(fiz), ^^"+1 = 1 + Q,2 +^2 'f" - 1 + ^2 + ^2 2i?2 '^"-' 

^ { " J^„(^/Vl + «^ + i8^ /?) + . ^ j„(«/Vl+a^ + ^\ «)) • (B8) 

Here the7„(-y, /?)'s used in (Bl) and (B8) are defined by 

jo = S'{yR), 

J2n+1 "= BJ2n ~ TlAJ2n-\ ~ «^n, 
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J2n = BJ2n-i " (/I " h)AJ2n-2, (B9) 

kn = -{n- m-'kn-u ko = iT-'l^yR-\ 

A = {l-y^)R-\ B = y\ 

In all the above equations, both a and /3 are assumed to be real and positive. However, these 
equations may also be used for complex values of a and jS if the proper formulas are used. These 
formulas may be found using analytic continuation. 
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